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We examine the violation of the majorization relations between the eigenvalues of the full and 
reduced density operators of entangled states of composite systems and its detection using gener¬ 
alized entropic forms based on arbitrary concave functions. It is shown that the violation of these 
relations may not always be detected by the conditional von Neumann and Tsallis entropies (for 
any q> Q). Families of smooth entropic forms which are always able to detect such violations are, 
however, provided. These features are then examined for particular sets of mixed states in a two- 
qudit system, which for d > 3 may exhibit different types of violation of the majorization relations. 

Comparison with the Peres criterion for separability is also shown. 
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Entanglement is now recognized as one of the most 
important features of quantum mechanics. It represents 
the possibility of composite quantum systems to exhibit 
correlations which cannot be accounted for classically. It 
plays therefore a fundamental role both in philosophical 
problems concerning the interpretation of quantum the¬ 
ory Q as well as in applications like quantum computa- 
tiom quantum cryptography and quantum teleportation 
[2-u| , where it is actually viewed as a resource. Formally, 
a quantum system composed of two subsystems A and B 
is said to be in a separable or unentangled state if the 
density operator of the whole system can be written as a 
convex combination of uncorrelated densities @ , 

p = '^LOaPA® Pb^ 0<a;Q<l, (1) 

OL 

where = 1 and p'^, are density operators 

for each subsystem. Otherwise, the system is said to 
be entangled or inseparable. When separable, the sys¬ 
tem fulhlls all Bell inequalities and satisfies other proper¬ 
ties characteristic of classical composite systems, making 
it unsuitable for the previous applications. Pure states 
(p = |IE')('I'|) are separable only for tensor product states 
(Ilk) = \(j)A)\4>B)), but in the case of mixed states it is in 
general much more difficult to determine whether a given 
density is separable or not. There are several simple nec¬ 
essary conditions for separability nn but they are not 
sufficient in general. 

Amon g th em, one is directly related with disorder and 
entropy [ij. Let 

Pa = Tib[p], pb = Tr a[ p], 

be the reduced density matrices of the subsystems, where 
Tr^.B denotes partial trace. It has been shown in [l^ 
that if p is separable, then p is more mixed (or disor¬ 
dered) than both pA and ps (the disorder criterion for 
separability). This property, which can be written as 
p < PA, p < Pb means in the present situation that 

i i 

Si = '^p^<S^ = '^pf, for f = l,...,dA, (2) 

i=i i=i 


(and similarly Si < Sf for i = 1 ,..., c^b) where pi, pf 
are, respectively, the eigenvalues of p and pA sorted in 
decreasing order and dA is the dimension of pA- The set 
of eigenvalues of p is then majorized by those of pA and 
Pb- For a pure state, Eqs. m can be fulfilled for A and 
B only ii Pi = pf = 1, implying that the state must 
be a tensor product state. However, for mixed states the 
disorder criterion is only a necessary condition for sep¬ 
arability. No sufficient condition for mixed states can 
actually be based on the knowledge of the eigenvalues of 
p and pA,B alone [3,[3. Nevertheless, it is the strongest 
spectral criterion and expresses a fundamental classical 
property of separable systems: a classical bipartite sys¬ 
tem is always more disordered than its subsystems, as 
the set of joint probabilities {pij} are always majorized 
by the set of marginal probabilities {pf = '^jPij}- 
The disorder criterion admits a simple entropic formu¬ 
lation [ 3 . Let us define the general entropic forms [13 


5/(p)=Tr/(p), 


(3) 


where / is an arbitrary smooth concave function satis¬ 
fying /(O) = /(I) = 0. These forms satisfy most basic 
properties of the conventional entropy except those re¬ 
lated with additivity [ 13 . In particular, if p ^ o', then 
Sf{p) > Sf{p') for any / of the previous form [ig. [l8l|. 
The disorder criterion implies then that in a separable 
state, the generalized entropy of the whole system is not 
less than those of the subsystems: 


-S'/(p) > Sf[pA), (4) 

and similarly, Sf{p) > Sf{pB)- In particular, the stan¬ 
dard von Neumann entropy 


‘S'(p) =-Trplnp, (5) 

corresponds to /(p) = —plnp, and Eq. ([4]) implies that 
the von Neumann conditional entropy 5'(p) — S{pa) [13 
is non-negative in any separable state |ll| . Similarly, the 
Tsallis entropy [l^ 


5,(p) = (l-Trp«)/( 9 -l), 


( 6 ) 
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corresponds to f{p) = (p — p'^) /(q — 1), which is concave 
for all q > 0, and Eq. ([2]) implies that the Tsallis condi¬ 
tional entropy [2^ [Sq{p) — Sq{pA)]/T'rp\ is non-negative 
in any separable state. For g —>• 1, Sq{p) —>■ S{p). 

For a given choice of /, Eq. (|4]) provides therefore a 
necessary test for separability of mixed states. In the 
von Neumann case, the criterion is actually rather weak, 
but in the Tsallis case the validity of Em (HJ for all g > 0 
provides a more stringent requirement 20l-l22l| , which has 
been shown to yield sufficient conditions for separability 
for some classes of states, like Werner states for n-qubits 
[ 2 ^ as well as any two-qubit state diagonal in the Bell 
basis 

Let us come now to the central point of this work. It 
is also true that if Sf{p) > Sf{p') for any f of the pre¬ 
vious form, then p < p' as will be explicitly shown 
below. Hence the validity of Eq. 0 for all f (and its 
partner for ps) provides an equivalent formulation of the 
disorder criterion. On the other hand, its validity for a 
particular / does not imply of course that p ^ pA- More¬ 
over, the validity of Eq. for all q > 0 in the Tsallis 
case does not imply p ^ pA- This may happen whenever 
the first violation of Eqs. 0 occurs for i > \. For suffi¬ 
ciently large q, Sq{p) Ri (1 — mip\)/{q — 1), where mi is 
the multiplicity of the largest eigenvalue pi of p. Hence, 
if Pi > Pi ^ Sq{p) < Sq{pA) for sufficiently large q. How¬ 
ever, if the first violation takes place for i > 1, Eq. © 
may remain valid for all g > 0 , as can be seen in simple 
examples. While in a two qubit system the only possible 
breakdown of Eqs. 0 is for i = 1 (as S 2 = 1), this is 
not the case in systems of larger dimensionality, where 
it is actually not rare to encounter entangled states in 
which the first violation of Eqs. 0 occurs for any i > 1, 
as shown in the example below. In such cases the Tsallis 
conditional entropy may not detect entanglement even 
when the disorder criterion does. The same will hold for 
any other set of functions f{p) whose maximum is ob¬ 
tained for p —^ I in the limit where it becomes sharp 

0 - 


We shall now provide a set of smooth concave functions 
which are always able to detect the breakdown of the in¬ 
equalities 0 , such that Sf[p) < Sf{pA) for some / of 
this set if Eqs. 0 are violated for some i > 1. This can 
be achieved with functions having a maximum at an op- 
timizable point a G ( 0 , 1 ), which can be made arbitrarily 
sharp by adjusting a second parameter. A possibility is 


fip) = 9t{p - a) - (1 - p)gt{-oi) - pgt{l - a ), (7) 

where a G [ 0 , 1 ] and gt{x) is a smooth concave function 
satisfying 


lim gt{x) = -\x\/2 . 


S 
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FIG. 1. The entropic functions 0 with g given by Eq. ®, 
for a = 1/4 and different values of t. fm = ce{l — a) denotes 
the maximum value in the t —>■ 00 limit. 


For any of these functions, /(O) = /(I) = 0, while for 
t —>■ 00 and 0 < a < I, 


f{p) 


p(I - a) , 
a{l-p ), 


p < a 
p > a 


(9) 


In this limit, 

Sf{p) = Tr[gt(p - a)] - {D - l)gti-a) - gt{l - a) 

a{na - 1 ) -I- 1 - '^Pij ( 10 ) 

Pi>a 

where D = dAds is the dimension of p and Uq, is the 
number of eigenvalues of p larger than a. It is now seen 
that if p -/ PAi Sf{p) < Sf{pA) for some a and t. Suppose 
that the inequalities 0 are first violated for i = j, i.e., 
Sj > which implies pj > p^. Choosing a = p^, we 
have Uaip) = j' > j, obtaining, for t —>■ 00 , 

Sfip) - SfipA) ^ pfif - j) - {Sp - Sf) 

< -[(Pf - pf)if - j) + Sj - ' 5/1 < 0 . 

Hence Sf{p) < Sf{pA) for sufficiently large t and a suf¬ 
ficiently close to pf. 

In the opposite limit t —>■ O'*", and assuming gt{x) = 
h{tx)/t, as in Eq. ([5]), f{p) —>■ t\h''{0)\{p — p^)/2, be¬ 
coming hence independent of a and proportional to the 
Tsallis case g = 2. On the other hand, for a = 0 or 1 
and t —>■ 00 , the lowest non-zero order of / depends on 
the higher order term rt(x) = gt(x) -I- |a;|/2 — gt(0). For 
instance, if a = 1 , f{p) = rt(p — 1 ) — (1 — p)rt(— 1 ), with 
rt{x) = —(2t)“^ ln[(l -I- in the case ([5]). 

Example. Let us consider a system of two qudits {d- 
dimensional subsystems) with d > S. We will first ex¬ 
amine a mixture of two antisymmetric states with a full 
random state, 


An example is (see Fig. 1) p = a;i|01 )(01 |-l-a; 2|02 )(02 \+ylA(SilB, (H) 

where \ij-) = {\ij) - \ji))/V2, y = {l-xi- Xa)/^^ and 
( 8 ) Ia,b = J2i=o K)(*l' X 2 = 0 and d = 2 , it reduces to 


gt{x) = — ( 2 t) ^ Incosh(te). 
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a Werner-Popescu state for two qubits [1, . Its eigen¬ 

values are obviously Xi-\- y, X 2 + y and y {{(P — 2 )-fold 
degenerate). Positivity of p imposes then the conditions 
xi + X 2 < 1 , Xi < {(P — l)xj -|- 1, i ^ j, which determine 
a triangle R with vertices at x = (xi,X 2 ) = (1, 0), (0,1), 
— (l,l)/((i^ — 2), shown in Fig. [2] for d = 3. Negative 
values of Xi are in principle feasible and represent a de¬ 
pletion of the antisymmetric states with respect to the 
fully mixed state. 

The lowest eigenvalue of the partial transpose of p is 

a{x) =y-\x\l2, ( 12 ) 

where |x| = {x\ -\- x\)^P. The equation aix) = 0 deter¬ 
mines an ellipse centered at x\= X 2 = —^l{(P — 8 ), with 
axes rotated 45 degrees with respect to the xi,X 2 axes, 
and cr(x) < 0 outside this ellipse. Hence, according to 
the Peres criterion Q , p will be entangled in the intersec¬ 
tion of this region with R (see Fig. [ 21 ). Moreover, it can 
be easily shown that p is indeed separable if a{x) > 0 . 
We may rewrite (fTTl) in the form 


n n 


p = '^Pi + y{iA ® Ib - y^^Qi) , 

i—1 i—1 

(13) 

Pi = Xi|0i“)(0i“| -1- yQi, 

(14) 


Q^ = (Zz|00)(00| + \0i){0i\ + K0)(f0| + \ii){u \, (15) 


where n = 2, qi > 0, qi + q 2 = 1. For xi,X 2 G R, pi is a 
non-negative operator for a two qubit system, where the 
Peres criterion is sufficient 0 , and will then be separable 
if and only if the lowest eigenvalue of its partial transpose 
is non-negative. This leads to the condition 

xl < 4qty^ ■ (16) 

If Eq. (fT6l) holds for i = 1,2, then |xp < 4y^, implying 
a{x) > 0. Moreover, choosing qi = xf/lxp, it is seen 
that any state satisfying <j{x) > 0 will fulfill Eqs. (HU) 
for i = 1 , 2 , implying separability of pi and hence of p. 
Thus, p is entangled if and only if a{x) < 0. Denoting 
with 7 the angle between x and v = ( 1 , 1 ), this implies 
entanglement for 

|x| > [^cos7-bdV2]“^ (17) 

The entanglement threshold for |x| is minimum at 7 = 0 
{xi = X 2 ), and decreases as d~‘^ as d increases (for d = 3 
the entangled region already covers 87% of R). It X 2 = 0 
Eq. (fTTl) implies entanglement just for 

xi > {1 + d^/2)-K (18) 

For d = 2 this yields |a;i| > 1/3, in agreement with the 
well-known result for a Werner-Popescu state 00 - 
As seen in Fig. [21 the disorder criterion detects entan¬ 
glement only in a sector of the entangled region (covering 
77% of it for d = 3), which can be classified according 
to the subset of inequalities ([2]) that are violated. Note 



FIG. 2. Region of allowed values of xi, X 2 in Eq. (Hill , for d = 
3. p is entangled in sectors a, b, c, d, and separable in e. The 
disorder criterion detects entanglement just in sectors o, 6 , c, 
where the inequalities m that are not fulfilled are indicated. 


that the complement of this sector is not convex. The 
reduced density matrix of any of the subsystems is 

PA = ^4^|0)(0| + ^|1)(1| + ^|2)(2| +ydlA. 

Hence, for Xi > X 2 > 0, 5i > if Xi+y > (a;i-|-X2)/2-|- 
yd, i.e.. 


Xi > 


1 + X 2 ((5 — 1) 

sTi 


~ 2{d-l) ’ 


(19) 


while if Xi > 0 > * 2 , this occurs for Xi + y > xi/2 + yd, 
i.e., xi > {1 — X 2 )/{t-\-5). Besides, for Xi > X 2 , S 2 > S 2 
in R just for xi -I- *2 + 2?/ > Xi -|- X 2/2 -|- 2yd, i.e.. 


xi>l-X 2 {l + 5/2). (20) 


These regions partially overlap, giving rise, together with 
their symmetric partners for X 2 > xi, to the sectors a, b, c 
of Fig. 121 In c, only the second inequality is violated. For 
3 < i < d, Si < Sf V xi, X 2 S R, with St = Sf only along 
the border xi -I- X 2 = 1 (where 5^ = 1 for i>2). 

In particular, for X 2 = 0 the disorder criterion predicts 
entanglement just for 

xi>(l + (5)-i, (21) 


which provides a bound that is strictly larger than Eq. 
m except for d = 2, and decreases only as d ^ as 
d increases. The same occurs for xi = X 2 (7 = 0), 
where Eq. (I2Q1) predicts entanglement just for |x| > 

[V 2 + 6 /{ 2 V 2 )]-\ 

For any entropic function /, the region where Eq. (HD 
does not hold, i.e., where 

Sfip) = Sfip) - SfipA), 


( 22 ) 
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FIG. 3. Top: Curves where Sf{p) = 0 according to: the 
Tsallis entropy for the indicated values of q (solid lines), and 
the entropy ([7]) (dashed line) with g given by Eq. ([8)l and 
a = 0.28, t = 30. The curve for g = 1 corresponds to the von 
Neumann entropy. Dotted lines indicate the boundaries of the 
sectors of Fig. [5] Bottom: Union of sectors where Sf{p) < 0 
according to: the von Neumann entropy (q = 1), the whole set 
of Tsallis entropies (all q > 0) and the whole set of entropies 
llTlIgll (t > 0, 0 < a < 1). 


is negative, will be necessarily contained in sectors a, b, c. 
As seen in Fig. [3] for d = 3, this region is rather small 
for the von Neumann entropy (<? = 1), but in the Tsallis 
case it increases as q increases, covering sectors a and b 
for q —>■ c», although it never covers sector c. Even the 
region where S^{p) < 0 for some g > 0 , depicted in the 
bottom panel, leaves out an appreciable fraction of sector 
c. The border of this region in sector c corresponds to 
finite and varying values of q, which can be determined by 
the simultaneous conditions S^{p) = 0, dSf{p)/dq = 0. 
Note also that along the outer border {xi + X 2 = 1), 





0.2 0.4 a 


FIG. 4. The behavior of Sf{p) for xi = X 2 = x (sector c of 
Fig. [3]| at the indicated values of x, according to: the Tsallis 
entropy as a function of q (top), and the entropy © for the 
case ®, as a function of t (center) and a (bottom). 

Sfip) < 0 for any f since here p"^ p. As q —^ 0+, the 
region where S^{p) < 0 shrinks and approaches this line. 

The typical behavior of the Tsallis conditional entropy 
Sf{p) = Sf{p)/Tr[p‘ff\ in sector c is depicted in the top 
panel of Fig.|3]as a function of q. Here Sf{p) may become 
negative only in a finite interval of q values. For xi = 
X2 = X, this occurs just for x > 0.381 (at this point 
Sf{p) — 0 just for q = 3.575). In comparison, along 
this line the disorder criterion predicts entanglement for 
a: > (2 + 5/2)~^ = 0.32, while the true value given by 
the Peres criterion is a: > (2 + df k, 0.12 (in the 

von Neumann case, S^{p) < 0 just for x > 0.452). Note 
as well that S^{p) does not necessarily develop a local 
minimum as a function of q in this sector (for xi = X 2 = x 
this occurs just for x > 0.334). 

In the case of the functions ((71), S^{p) becomes nega¬ 
tive for sufficiently large t at any point of sectors a, b, c 
if a is appropriately chosen. Due to the symmetries of 
the present case, just two values of a are required: For 
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d = 3, in sectors a, 6 it is sufficient to choose any single 
a > 5/13 ~ 0.385, which is the value of pf at the tip of 
sector b (x = (4/13,0)), whereas the whole sector c can 
be covered choosing a = 7/25 = 0.28, which is the value 
of P 2 at the inner tip of this sector {xi = X 2 = 0.32). In 
the case ([ 8 ]), the curves where Sf{p) — 0 are, for large t 
and a > 1/2, similar to those of the Tsallis case for large 
q, whereas if a = 0.28, they approach for large t the in¬ 
ner border of sector c, as seen in the top panel of Fig. [S] 
The behavior of ^(p) = Sf{p)/\TT[gtip)]\ in sector c is 
depicted in Fig. 01 where it is seen that for sufficiently 
large t, it becomes negative in a finite interval of a values 
(which always includes a = 0.28), for any x > 0.32. 

Similar considerations hold for the states 

n 

p = '^Xi\0i~){0i~\+ylA<Si Ib , (23) 

where n < d—1 and y = (1 — Xi)fd'^. Its eigenvalues 
are Xi +y, i = 1,... ,n, and y ((d^ — n)-fold degenerate), 
and the allowed values of Xi are contained in a region R 
bounded by the conditions y > 0, Xi > —y. 

The lowest eigenvalue of the partial transpose is again 
given by Eq. (IT^ . with \x\ = and it can be 

seen that p is entangled if and only if a{x) < 0 : writing 
p in the form (fT3l) . with qi > 0, Ym=i di = 1) state 
satisfying cr(a;) > 0 is seen to fulfill Eqs. (fTCl) for i = 
1 ,... ,n if qi = xf/|a:p, which implies separability of pi 
and hence, of p. The entanglement condition ij{x) < 0 
can be cast as 

\x\ > (-v/ncosy -I- d^I2 )~^, (24) 

where 7 is the angle between x and v = (!,...,!), and 
corresponds to the region outside an n-dimensional el¬ 
lipsoid. For I 7 I < 7 r /2 and fixed d, the entanglement 
threshold for \x\ decreases as n increases. 

For n > 5 (and d > n -I- 1) a novel feature arises: en¬ 
tanglement may also occur for Xi < 0 V i, i.e., in the 
vicinity of 7 = tt, since the vertex of R in this direction 
{xvi = — (d^ — n)“^ for i = 1 ,..., n) lies outside the el¬ 
lipsoid: |xt,| = y/nj{d‘^ — n) > {d? 12 — y/n)~^ for n > 4. 
Thus, for n > 5 entanglement can also occur by deple¬ 
tion of antisymmetric states, although the corresponding 
region is very small (see Fig. [5]for the case n = 5, d = 6 ). 
Nevertheless, the ellipsoid < 7 ( 0 ;) > 0 is never fully con¬ 
tained in R. Points of the boundary of R like Xj = 0 if 
j i, Xi = —l/(d^ — 1), fulhll a{x) > 0 V d > 2. 

The disorder criterion ensures again entanglement in a 
much smaller region. Moreover, the first violation of the 
inequalities ([ 2 ]) may now occur for any i between 1 and 
n. The reduced density of system A is 

1 " 

PA = **[l*)(*l + |0)(0|] + ydiA ■ (25) 

For Xi > X 2 > ... > x„ > 0 and 1 < i < n, Si > for 

1 "■ 

Xi > +'iyid- 1 ). 

i=i 



FIG. 5. Section of allowed values of x, X 5 in the state (1231) for 
d — 6 , n = 5 and Xi = x for i = 1,..., 4. p is entangled in all 
shaded sectors. The disorder criterion detects entanglement 
in sectors a, b, c, d, e, where the inequalities © that are not 
fulfilled are indicated. Also depicted are the curves where 
^fip) = 0 in the von Neumann case (solid line) and the 
Tsallis case g = 2 (dashed line), such that Sf{p) < 0 above 
these curves. The dotted line indicates the boundary of the 
region covered by all g > 0 in the Tsallis case (see text). 

According to the values of the Xi, we may have violation 
of the inequalities © for any set of m distinct indices 
jj, with m < n, ij < n. In particular, if Xi = x > 0 
for i = 1,... ,n (7 = 0 ), only the inequality can be 
violated: > 5;/ for 

X > nf (n^ -h (5), (26) 

i.e., |a:| > [\/n + This bound is larger than 

that given by Eq. (l24ll for 7 = 0, for any n > 1, d > 3. 
On the other hand, the disorder criterion is not able to 
predict entanglement for x < 0 (7 = tt) for any d, n. 

Eor instance, Eig. [5] depicts, for n = 5 and d = 6 , the 
intersection of R with the plane Xi = x for i = 1,..., 4. 
The separable sector covers only a small fraction of this 
section (~ 2 , 6 %) but the disorder criterion covers only 
~ 51% of the entangled region. Moreover, the first vio¬ 
lation of the inequalities ([2|) may occur for f = 1, 4 or 5, 
the last two in sectors c, d, e, although the largest region 
corresponds to z = 1 (sectors a, b, accounting for 40% of 
the entangled region). 

The region where S^{p) < 0 is small in the von Neu¬ 
mann case, but increases as q increases in the Tsallis 
case, covering sectors a and & for g —>■ 00 , although the 
covering of sectors c, d, e is very poor. The inner border 
of the region where S^{p) < 0 for some g > 0 lies very 
close to the outer boundary (X]r=i = 1 ) these sec¬ 
tors, being indistinguishable from it in the scale of the 
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figure. For instance, for X5 = x (sector c), the bound¬ 
ary lies at X = 0.2 and S^{p) < 0 for some q just for 
X > 0.19997, while for X 5 = 0 (sector e), the boundary is 
at X = 0.25 and Sf{p) < 0 only for x > 0.2492. In com¬ 
parison, along these lines the disorder criterion predicts 
entanglement for x > 0.175 and x > 0.204 respectively, 
as given by Eq. ®. For the entropy 0, it is necessary 
to employ in this case an interval of a values in order to 
have Sf (p) < 0 for large t in all points of regions c, d, e 
(a S [0.069,0.108] in c and a. G [0.114,0.134] in e), al¬ 
though for the whole sectors a and b it is sufficient to take 
a single a > 1/4. It should be mentioned that all previous 
equations and results remain valid if the antisymmetric 
states in (HI) or (l23l) are replaced by symmetric states 
= (l*j) + |j*))/'/2- 


In conclusion, we have provided an example of a set 
of smooth entropic forms which are always able to de¬ 
tect the breakdown of the majorization relations ([2|) in 
entangled states, through the sign of the entropic differ¬ 
ence (I 22 I) . We have also shown that when the first viola¬ 
tion occurs for i > 1, the Tsallis conditional entropy may 
not detect such violations for any q, being hence weaker 
in these cases than the disorder criterion. The present 
entropic forms are therefore necessary for detecting and 
measuring entanglement by entropic means in such sit¬ 
uations. These issues have been illustrated in detail for 
particular states of a two qudit system, which show that 
the first violation of Eqs. may in fact occur for any 
i < d — 1. 
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